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We establish a connection between non-Markovian memory effects and thermodynamical quantities such as 
work. We show how memory effects can be interpreted as revivals of work that can be extracted from a quantum 
system. We prove that non-Markovianity may allow an increase in the extractable work even when the entropy 
of the system is increasing. Our results have important implications both in quantum thermodynamics and in 
quantum information theory. In the former context they pave the way to the understanding of concepts like 
work in a non-Markovian open system scenario. In the latter context they lead to interesting consequences for 
quantum state merging protocols in presence of noise. 


INTRODUCTION 

The connection between thermodynamics and information 
theory, expressed by Landauer’s principle, is a milestone of 
the physics of the last century. According to this principle, the 
erasure of information stored in a system requires an amount 
of work proportional to the entropy of the system. The nat¬ 
ural framework to discuss thermodynamics at the quantum 
level is the theory of open quantum systems. A number of 
recent results have shown that memory effects arising from 
strong system-environment correlations may lead to informa¬ 
tion back-flow, hence prolonging the life of quantum prop¬ 
erties. Open systems exhibiting such behaviour are known 
as non-Markovian. The relation between non-Markovianity 
and quantum thermodynamics has been until now largely un¬ 
explored. Flere we establish this missing link by means of 
Landauer’s principle. We show that memory effects control 
the amount of work that one can extract from an open quan¬ 
tum system. Hence, the work extraction can be optimised via 
reservoir engineering techniques. 

Landauer’s erasure principle states that, in order to perform 
irreversible operations on a system, such as erasing a bit of in¬ 
formation, a certain amount of work must be performed. This 
is known as the work of erasure IH. The very same princi¬ 
ple holds in the quantum regime canning far-reaching conse¬ 
quences Hi. As noticed in Ref. HI, indeed, in a quantum 
setting the presence of a quantum memory correlated to the 
system leads to the possibility of extracting work while eras¬ 
ing information on the system. This possibility crucially relies 
on the persistence of quantum correlations between the system 
and the quantum memory. One would therefore expect that, 
due to the usually deleterious effect of the environment acting 
on either the system or the memory, this uniquely quantum 
feature should only survive for a very short time. In other 
words we expect that, after the initial preparation of a quan¬ 
tum correlated state between system and memory, the ability 
to extract work by erasure will characterise only a short initial 
transient. 


One of the main results of this Article is to prove that this 
needs not be true for non-Markovian systems Od UlTIl . We will 
show how memory effects are directly connected to the tempo¬ 
ral behaviour of the work of erasure and are therefore instru¬ 
mental to the persistence of work extraction. This is impor¬ 
tant also in the ongoing discussion on the usefulness of non- 
Markovianity for quantum technologies ifoUl^ . since our re¬ 
sults prove for the first time that memory effects are useful for 
work extraction purposes. 

LANDAUER PRINCIPLE IN THE OPEN SYSTEM 
SCENARIO. 

The extension of Landauer’s principle to quantum systems, 
formulated by Lubkin, states that the work performed on a 
quantum system S to erase classical information stored in it, 
is connected to the system’s entropy via the relation 

W{S) ^ H{S)kTln2, (1) 

where H(S ) is the von Neumann entropy quantifying our lack 
of knowledge on the system, k is Boltzmann constant, and T 
is the temperature of the reservoir used for the optimal erasure 
process. 

When one tries to ascertain properties of the system, how¬ 
ever, one is naturally led to introduce an observer O which 
performs measurements on 5. In this framework the “amount 
of information”stored in a system is observer-dependent. One 
should therefore replace the von Neumann entropy H{S) in 
Eq. ([T]i with the conditional entropy H{S\0). The fully quan¬ 
tum scenario is described by a situation in which the observer 
possesses a quantum memory Q. In Ref. |[3l it is shown that, 
in the thermodynamic limit, the work of erasure is given by 

W{S\Q)^H{S\Q)kT\n2, (2) 

and the extractable work, that is the work that one can extract 
from an «-qubit system, is 

We,^[n-H{S\Q)]kT\n2, (3) 
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FIG. 1. (Colors online) Illustration of the two physical scenarios 
considered in the Article: a) the open quantum memory and b) the 
open quantum system. 


where H(S\Q) - H(S Q) - H{Q) is the conditional entropy, 
i.e. the entropy of the system conditioned on the memory, 
with H(S Q) and H(Q) the von Neumann entropies of the com¬ 
bined system-memory and of the memory, respectively. If the 
system and the memory are prepared in a quantum correlated 
state, such as maximally entangled state, the conditional en¬ 
tropy may be negative. When this happens one can extract 
work while erasing information on the system. This is a gen¬ 
uinely quantum phenomenon with no classical counterpart. 

We will now extend this analysis to the case in which ei¬ 
ther the system or the memory is subjected to the action of the 
environment. This will give us a basic understanding of the ro¬ 
bustness of the quantum-enhanced work extraction to realistic 
and unavoidable sources of noise. The two situations are pic- 
torially illustrated in Fig. 1. Both the work of erasure and the 
extractable work are in these scenarios time-dependent. One 
would expect that, due to the presence of the environment, the 
work of erasure would increase in time while the extractable 
work would decrease. Let us indicate with 


AW,,(fi,f2) = W„(f2)-W„(fi) 

= [H(S\Q)r,-H{S\Q\]kTln2, (4) 

with f 2 > fi, the variation in time of the extractable work, 
and with AW(fi, 12 ) = -^WexihJi) the change in the work of 
erasure. 

It is worth stressing that the processes of work extraction or 
erasure IBfil are not directly considered in our scheme and 
are assumed to take place at a given time f, after the system 
has interacted with the environment. The optimal process to 
extract the amount of work given by Eqs. (l2]i-([3]) is described 
in |13l ■ Our main goal is to understand how the environment 
affects, as time evolves, the ability to extract work from the 
system. By using reservoir engineering techniques 123-221 
one may then be able to minimize the detrimental effects of the 
environment. As we will see, memory effects will be the key 
ingredient to accomplish this goal. We begin by considering 
the case of an open quantum memory. 

Open quantum memory: The action of the environment on 
the quantum memory is described in terms of a dynamical 
map <!),, that is a f-parametrized family of completely posi¬ 
tive and trace preserving (CPTP) maps. We indicate the con¬ 
ditional entropy in Eq. (|4]i with H{S\Q,), with Qt - Of(2), 
to emphasize that only Q is evolving due to the coupling to 


the environment. This quantity is just the negative of the co¬ 
herent information, i.e., I{S)Qt) — -H{S\Qt), measuring the 
correlations between the system and the memory ll23[l . Hence, 


AW„(fi, f2) = [KS )Q,,)- I{S)Q„ )] kT In2. (5) 

Eor divisible dynamical maps, namely when d), = with 

s < t and Of j CPTP, the data processing inequality implies, 
that 


KS)Q,,) < I{S)Q„), forfi<f2. (6) 

In this case, the extractable work monotonically decreases in 
time and the work of erasure monotonically increases, as one 
would intuitively expect. 

Open quantum system: Let us now investigate what hap¬ 
pens when the environment acts on the system rather than on 
the memory. The extractable work is now given by 

Wexit)^[n-H{S,\Q)]kTln2 

= [n-H(S,) + I{S,:Q)]kT\n2, (7) 

where H{St\Q) indicates the conditional entropy when the 
system is coupled to the environment, i.e., St - d),(5'), and 
I (St ■ Q) - H(St) + H(Q) - H(StQ) is the mutual informa¬ 
tion at time t, quantifying the amount of information that S t 
and Q share «2M . 

The change in the extractable work can be recast in the form 

A W„(f 1 , f 2 ) = [-AH(S t) + AI(S ,:Q)]kTln2, (8) 

where A//(5,) = H(S„)-H(S „) and AI(S, : Q) = I(S„ : Q)- 
I(Sti ■ 2)’ t 2 >t\. Once again, for divisible dynamics, 

I(S„-.Q)<I{Stt-Q), forfi<f2, (9) 

which implies that an increase in extractable work requires a 
decrease in entropy of the system. This meets our intuition, 
since it amounts to saying that the more we know about the 
system, namely the closer it is to a pure state, the more work 
we can extract from it. Equivalently, the closer the system is to 
a maximally mixed state, the smaller is the extractable work. 
As we will see in the following, however, non-Markovian 
memory effects may defy this intuition. 

We note in passing that for unital dynamical maps the 
change in the entropy of the system vanishes when the initial 
system-memory (2 qubit system) is prepared in a maximally 
entangled state. In this case we obtain exactly the same result 
of the open quantum memory scenario, that is, the variation in 
the extractable work is proportional to the variation in the co¬ 
herent information, see Eq. ©, since AI{St : Q) = AI{S,)Q). 


THERMODYNAMIC MEANING OE NON-MARKOVIANITY 

The most general form of open system evolution is de¬ 
scribed by dynamical maps which may violate the divisibility 
property while still being CPTP. When this happens, the co¬ 
herent information may temporarily increase for certain time 
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FIG. 2. (Colors online) Time evolution of the extractable work 
WexIkT ln2 for one qubit system. The blue solid line corresponds 
to a Pauli channel model with decay rates ji(t) = 72(0 = ^/2, 
and 73(0 = tu tan(tuf)/2, with T = 0.1 and tu = 2 (a.u.). The red 
dashed line corresponds to a Pauli channel model with decay rates 
7i(0 = 72(0 = ^/2, and 73 (f) = -wtanh(tuf)/2, with A = 1 and 
a> = 0.5 (a.u.). 


intervals. This, in turn, generally leads to a partial recovery of 
extractable work or, equivalently, to a decrease in the work of 
erasure. 

Recent research on open quantum systems has high¬ 
lighted the fact that memory effects can be associated to 
different manifestations of non-Markovianity which, in turn, 
are revealed by different physical quantities. Several non- 
Markovianity measures capturing this plethora of phenom¬ 
ena have been introduced and used to describe the informa¬ 
tion back-flow and the recoherence iini. A connection be¬ 
tween memory effects and thermodynamical quantities, how¬ 
ever, was until now missing. Equations (0 and ([ 8 ]l bridge this 
gap by showing how memory effects are indeed directly re¬ 
lated to work. 

Until now our results are model-independent, hence fully 
general, but formal. The key relevant question is: are there 
physical models of realistic environments in which the ther¬ 
modynamic power of non-Markovianity can be observed and 
tested? We will positively answer to this question by consid¬ 
ering three exemplary cases, namely a qubit coupled to a clas¬ 
sical environment (Pauli channel), a spin environment (Ising 
model) and a bosonic environment (amplitude damping chan¬ 
nel). 

In our first example we study the time evolution of the ex¬ 
tractable work for time-dependent Pauli channels. Pauli chan¬ 
nels describe depolarising noise acting independently along 
the X, y, and z directions with time-dependent rates 71 (f), 72 (f) 
and 73 (f), respectively. In this case, due to the unitality of 
Pauli channels, the behaviour of extractable work in time is 
described by Eq. (|5]) both in the open quantum memory and 
in the open quantum system scenario. In Eig. 2 we compare 
two different open dynamics. The corresponding behaviours 
of the extractable work are given by the solid blue line and 
the dashed red line in Eig. 2. In both cases the dynamical 


map is non-divisible, but only in one case oscillations of ex¬ 
tractable work do take place. When these oscillations occur, 
the uniquely quantum feature of negative conditional entropy 
may show revivals even after becoming positive due to noise 
induced by the environment. Eor detailed calculations of the 
time evolution of the extractable work we refer to the Supple¬ 
mentary Material. 

This example illustrates that, when defined in terms of co¬ 
herent (or mutual) information |@], the physical interpretation 
of non-Markovianity as information back-flow carries along a 
powerful connection to revivals of the extractable work. We 
can therefore interpret memory effects in terms of a quan¬ 
tity with a clear physical meaning, that is work, rather than 
the somewhat elusive concept of information. Equation (|5]), 
indeed, straightforwardly illustrates that memory effects are 
manifested as time revivals of work. 


NON-MARKOVIAN ENHANCEMENT OF WORK 
EXTRACTION BY ERASURE 

As the second example, we consider a physical system con¬ 
sisting of a central spin coupled to a spin chain. More pre¬ 
cisely we focus on the Ising model in a transverse field 1241] . 
Details on the Hamiltonian and the exact solution for the spin 
system are given in the Supplementary Material. This systems 
exhibits a quantum phase transition between ferromagnetic 
state and paramagnetic state when the parameter A*, measur¬ 
ing the ratio between the strength of the transverse field and 
the coupling between the environmental spins equals unity. In 
Ref. i25n it has been shown that the central spin dynamics is 
non-Markovian for any value of A*, except at phase transition. 
This is therefore an ideal system to study the modification of 
the dynamics induced by non-Markovianity as we can control 
the non-Markovian character via the parameter A*. We notice 
that, also in this case, the dynamical map is unital. 

Once again, due to unitality, the extractable work by erasure 
is just proportional to the coherent information both in the 
open system and the open memory scenarios. Fig. 3 shows the 
behaviour of extractable work for three exemplary values of 
A*. For d* = 1 (black solid line) the extractable work quickly 
decreases to unity signalling that the erasure work goes to zero 
after a short transient time. In this case the spin environment 
is at quantum phase transition and the dynamics of the central 
spin is Markovian. 

For A* =0.1 (green short-dashed line), after decreasing to 
unity, the extractable work presents oscillations with maxi¬ 
mum amplitude damped in time. These oscillations, corre¬ 
sponding to revivals of the coherent information, characterize 
the dynamics when the environment is in its ferromagnetic 
ground state. Finally, for A* = 1.9 (orange long-dashed line) 
oscillations in the extractable work are also present while this 
quantity remains larger than unity at all the times, meaning 
that work extraction from the erasure is always possible. This 
corresponds to the case, in which the environment is in a para¬ 
magnetic state. This example illustrates how non-Markovian 
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FIG. 3. (Colors online) Time evolution of the extractable work of 
one qubit for an Ising model in a transverse field. We compared 
the dynamics for N = 4000 spins and three values of renormalized 
transverse field 4: 4* = 1 (black solid line), 4* = 0.1 (green short- 
dashed line) and 4* = 1.9 (orange long-dashed line). 


FIG. 4. (Colors online) Time evolution of the mutual information 
(blue long-dashed line), system entropy (red short-dashed line), and 
extractable work (solid black line), for the amplitude damping pho¬ 
tonics band gap model, with detuning from the band gap edge fre¬ 
quency d = -1 (a.u.). 


memory effects allow to retain or regain extractable work, 
leading to a thermodynamic enhancement compared to the 
Markovian case. 


THE POWER OF QUANTUM CORRELATIONS 

The third example highlights the interesting situation in 
which memory effects are responsible for a counter-intuitive 
behaviour of the extractable work. In the open quantum sys¬ 
tem scenario the change in the extractable work is given by 
Eq. (|8j, namely it is the difference between the change in mu¬ 
tual information and the change in entropy. When the open 
system dynamics is Markovian, the mutual information mono- 
tonically decreases. Hence, in this case an increase of the 
extractable work can only occur when the entropy of the sys¬ 
tem decreases. Non-Markovian dynamics, however, allows 
for revivals of the extractable work even when the entropy of 
the system is constant or increasing. Let us discuss these two 
cases separately. 

We begin by considering non-Markovian dephasing along 
an arbitrary direction, i.e., a Pauli channel with two null de¬ 
phasing rates and the third one time-dependent. For a generic 
initial state of the system one can easily prove that the entropy 
change is either null or decreasing whenever the mutual in¬ 
formation increases. Therefore the best scenario for memory- 
induced revivals of extractable work corresponds to the case 
in which the change in entropy is zero. This is the situation 
occurring, e.g., in the Ising model in a transverse field, see 
Fig. 3. 

Let us now focus on the more counter-intuitive situation in 
which the increase in the coherent information is greater than 
the increase in the system’s entropy, that is, we get an over¬ 
all increase in the work extraction even though our knowl¬ 
edge of the system decreases. Since we assume that the initial 
system-memory state is maximally entangled, this situation 


may occur only for non-unital open system dynamics. Fig. 4 
illustrates the dynamics of the extractable work, as defined in 
Eq. (I8]l, the mutual information and the system entropy for an 
amplitude damping photonic band gap model. Once again all 
quantities can be calculated analytically; see the Supplemen¬ 
tary Material. 

In the figure, the green shaded regions highlight time in¬ 
tervals in which memory effects are so dominant to induce a 
partial recovery of the extractable work even if the system en¬ 
tropy increases. This apparent contradiction is in fact resolved 
when one thinks that, in this setting, what counts is the amount 
of information on the system possessed by an observer with a 
quantum memory. Therefore, it is clear that an increase in 
the system-memory correlations, indicated by the behaviour 
of the mutual information, plays a key role in determining the 
extractable work. 


APPLICATION TO QUANTUM INFORMATION THEORY 

We conclude by noticing that our results carry along inter¬ 
esting consequences when applied to one of the most impor¬ 
tant protocols in quantum information theory, namely quan¬ 
tum state merging. In quantum state merging two parties, Al¬ 
ice and Bob, share a common state pab, and Bob has to repro¬ 
duce this state having unlimited access to local operations and 
classical communication. The quantum cost of the protocol, 
i.e., the minimal amount of Bell pairs that Alice and Bob need 
to share in order to complete the task, is given by the condi¬ 
tional entropy H{A\B) 12611 . We see that the crucial quantity 
here is the same as in the work extraction protocol. Hence, 
one can think, analogously to what we did in the work extrac¬ 
tion protocol, in terms of the degradation of the quantum cost 
when either Alice or Bob are subjected to an environment. 

The same lines of thought that we used above to illus¬ 
trate the thermodynamic power of non-Markovianity lead us 
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to conclude that reservoir engineering can be used to induce 
memory effects hence lowering the quantum cost of the state 
merging protocol in presence of unavoidable noise. Moreover, 
it is known that quantum state merging is a primitive, i.e. it 
can generate all the other building blocks of quantum Shan¬ 
non theory 1271] . Hence reservoir memory effects are poten¬ 
tially beneficial for any quantum information protocol as they 
can improve their efficiency. 


CONCLUSIONS 


The interpretation of non-Markovian memory effects in 
terms of revivals of either the erasure work or the extractable 
work unveils a fundamental link between the theory of open 
quantum systems and quantum thermodynamics. Such a con¬ 
nection also gives a powerful answer to the question of the 
potential usefulness of non-Markovianity, when we think in 
terms of reservoir engineering. One of the current challenges 
of open quantum system theory is the development of a re¬ 
source theory of non-Markovianity. It is worth underlining 
that both in quantum thermodynamics and in quantum infor¬ 
mation theory resource theories have been successfully devel¬ 
oped 12811 . Our approach may therefore pave the way to estab¬ 
lishing, in the most general framework, if and in which way 
memory effects can be considered as a resource. 
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Supplementary Material 


In this Supplementary Material we present in detail the 
mathematical description of the physical models discussed in 
the main Article. In all the considered examples we take as 
an initial state of the system S and quantum memory Q the 
maximally entangled two qubit state, \SQ} - ^ (|10) - |01)). 
The motivation for such choice comes form the fact that the 
maximally entangled states are optimal states from the point 
of view of the work by erasure protocol |[ll]. 

In the Article we investigate two scenarios, when either the 
qubit of quantum memory or the qubit of system is subjected 
to environmental noise, hence undergoing a non-unitary evo¬ 
lution. In general, these two situations have different implica¬ 
tions on the behaviour in time of the extractable work (see the 
Eqs. (5) and (8) of the main Article, respectively). However, 
if the dynamical map is unital and the reduced state of the sys¬ 
tem is maximally mixed, as it is in fact the case in all studied 
examples, then the two scenarios are equivalent and provide 
a direct link between the revivals of the extractable work and 
the revivals of coherent information of the dynamical map. 

Two of the analysed physical models, Pauli channels and 
Ising model, are unital maps, which allows us to clarify the 
relation between divisibility, revivals of the coherent informa¬ 
tion and the work extraction. The last example, photonic band 
gap amplitude damping, is a non-unital map. Hence, we can 
study the more complex situation described in the Article’s 
Eq. (8), where also the system entropy is evolving, and wit¬ 
ness the surprising behaviour of the extractable work in time. 

Pauli Channel 

As the first example of a non-Markovian qubit channel, we 
consider a Pauli channel with dynamical map of the form 

3 

= X CTaPO-a, (1) 

ff =0 

where 0 < pa < 1 satisfy YiaPa - era are the Pauli 

matrices (cro = 1 ). A local generator L, corresponding to this 
channel is defined via 

3 

Uip) = ^ 7/(0 {o-ipo-i - p), (2) 

1=1 

where the relation between probability distribution Pa{t) and 
local decoherence rates ykii) is given by ||2l 

Po(0 = ^[1 + di(f) + A 2 {t) + / 13 (f)] 

Pl(f) = -[1 -t- /li(f) - A2(t) + /13(f)] 

PliO - -[1 - /ll(f) + A2(t) - /13(f)] 

P3(f) = ^[1 - di(f) - A2(t) + /13(f)], (3) 


and the time-dependent eigenvalues Ak(t) read 

/Ii(f) = exp(- 2 [r 2 (f) + r 3 (f)]), 

/ 12 (f) = exp(-2[ri(f) + r 3 (f)]), (4) 

/ 13 (f) = exp(-2[ri(f) + r 2 (f)]), 

with r^(f) = 7 k(T)dT. A Pauli channel is divisible if and 
only if YkCO satisfy the following conditions 

ri(f) > 0, 72(0 > 0, 73 (f) > 0, (5) 

for all f > 0. 

In a scenario of evolving quantum memory the extractable 
work is then given by 

W,At) = [l+HS}Qr)]kT In 2. (6) 

In this Pauli channel example, the coherent information 
I(S)Qt) for the initial state \SQ) - ^ (|10) - |01)) is easily 
analytically calculated: 

I(S >20 = 1 + X Iog2(p<.(0). (7) 

a 

In the Article we analyse, and depict in Pig. 2, the time evo¬ 
lution of the extractable work corresponding to two cases of 
Pauli channels with different values of the 7-rates. Por the first 
example we choose 71(f) - j - 72(f) and 73(f) = ^ tan(wf) 
that clearly violate the divisibility condition of Eq. © and 
show revivals in the extractable work (see the solid blue line 
in Pig. 2 of the main article). 

Por comparison, as the second example we choose the case 
analysed in ||2,[1], 71(f) - j - 72(f) and 73(f) = tanh(wf). 
A Pauli channel with such 7-rates also leads to legitimate non- 
divisible evolution whenever 0 < w < /I, however in this 
case no revivals of the extractable work are present (see the 
red dashed line in Pig. 2 ). This is to show that sheer non¬ 
divisibility of the evolving channel does not always imply in¬ 
creasing work gain. Interestingly, the second example cor¬ 
responding to 73(f) = -7 tanh(cuf) gives rise to so called P- 
divisible evolution |0] which is characterized by conditions 

7/(f)+ 73 (f) >0, ( 8 ) 

for i + j, i, j - 1,2, 3 and all f > 0. As was advocated in fl, 
a quantum evolution which violates Eq. ©, but satisfies Eq. 
®, displays only weak non-Markovianity. These examples 
show that essentially non-Markovian evolution, i.e., the one 
that violates both Eq. © and Eq. ©, may lead to increase in 
the extractable work, whereas weakly non-Markovian satisfies 
AlT„(fi,f2)<0. 

Ising Model 

As a second example we consider a qubit (Q or S) trans¬ 
versely coupled to an Ising spin chain in a transverse field ||3]. 
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Such interaction results in a pure dephasing dynamics of the 
qubit described by the following master equation |@] 

L,{p) = + y{t){a--j)a-^ - p), (9) 


with A(f) being the time-dependent Lamb shift and y{t) a 
time-dependent decay rate further related to the Loschmidt 
echo L{t) via relation 


r(0 = - 


m 

AL{ty 


( 10 ) 


Notice that this corresponds to a special case of a Pauli chan¬ 
nel with yi{t) = yiit) - 0 and y 3 (t) - y(t). Using the equa¬ 
tions 0 and 0 it is straightforward to show that 

/(5)a) = 1 + ^(1- Vi^log2(^(l- Vi«)) 

+i(i+Vi^iog2[^(i+(11) 

which can be inserted to Eq. 0 to obtain the extractable 
work. For this model the Loschmidt echo depends on two rel¬ 
evant parameters, A describing the strength of the transverse 
field and 6 the qubit-environment coupling strength. The ana¬ 
lytical expression for the Loschmidt echo is given by 

L{A, 6, t) = ni>o [l - sin^(2ySi) sin^(4o], (12) 


where are the Bogoliubov angles and the single quasi¬ 
particle excitation energies of the system with the qubit in the 
state |1) (for more details see Refs. Hi). 

In Fig. 3 of the Article we present three curves describing 
the time evolution of the extractable work for d = 1.8 (orange 
long-dashed line), A - 0.0 (green short-dashed line) and A - 
0.9 (black solid line). The remaining parameters are fixed for 
all the plots as 5 = 0.1, 7 = 1 and the number of spins N - 
4000. 


Photonic Band Gap 


The third physical model that we analyse is the amplitude 
damping channel described by the following exact local gen¬ 
erator: 

is{t) r 1 t 

Liip) =-^[cr+cr-.p] + y{t)[a-^po-+ - -{cr+cr_,p)j, (13) 


with s(t) = -2Im[G(f)/G(f)] and 7 (f) = -2 Re [G(f)/G(0] 
the time-dependent Lamb shift and decay rate, respectively. In 
Eq. (IT 3 I 1 cr+ and cr_ stand for raising and lowering operators, 
respectively, and the complex function G(f) satisfies |G(f)| < 
1 for all t > 0, and directly depends on the features of the 
reservoir. 

The dynamical map of this evolution can be described in 
terms of the Kraus representation <[),(p) = Ki(t)pKi(ty, 
where 


Kiit) = 


1 0 
0 G(f) 


,A2(f) = 


0 Vl - IG(0P 

0 0 


(14) 


We consider the model described in Ref. 101, in which the 
reservoir is an ideal photonic crystal with a frequency band 
gap. In such a case the specific form of G(f) is as follows: 

G(f) = 2vixie^-‘i'+'* + V 2 (X 2 +y 2 )e^"^'^'* 

- ^ ajyjll - (I)( (15) 

2=1 


where Ch is an error function. One can see that G(f) de¬ 
pends on two relevant parameters, the detuning 6 = cjq - a>e 
from the band gap edge frequency cOg, which we set to 5 = 
- 1 , and the characteristic frequency parameter fi defined as 
16neotic^ with eo the Coulomb constant and d 
the atomic dipole moment. In addition: 


xi = (A+ + 

1/2t1/3 



1 1 

2 “ 2 

4 6^ 



U 

{Xj - Xi)(Xj - Xk) 

yj = 

/ V -2 

(;■= 1,2,3) 


(16) 
(17) 

(j + i + k;j,i,k= 1,2,3), (18) 
(19) 


Notice that this dynamical map is not unital, so that one 
does not have the symmetrical situation with respect to the 
system and the memory dynamics, that was present in the 
previous models. Hence, it is a good candidate for investi¬ 
gating the counter-intuitive behaviour of the extractable work 
described in the Eq. ( 8 ) of the Article. We therefore focus 
on the case where the system is subjected to the amplitude 
damping channel described above. As in previous examples, 
the initial state is Q) - "^(110) ~ |01))- The entropic quan¬ 
tities entering in the definition of the extractable work can be 
calculated analytically 

H(S,) = (2 - |G(r)|2)log2 ((2 - |G(f)l')j 

-i|G(f)l'log2(i|G(f)pj, (20) 

H(Q) = 1, ( 21 ) 

H(S,Q) = - 5 ( 1 - |G(0l')log2 (1 - |G(f)l")j 

(1 + |G( 0 l")log 2 (1 + |G(f)l')j, ( 22 ) 


and further used to calculate the mutual information I(S , : 
Q) - H(St) + H{Q) — H(StQ). Using these we get the analyt¬ 
ical solution for the extractable work: 


= [1 - //(5,) + 7(5, : 0] krin2. (23) 

This example, where the system is subjected to the pho¬ 
tonic band gap amplitude damping, illustrates that the counter¬ 
intuitive situation in which one has revivals of extractable 
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work even though the entropy of the system increases, can 
occur in realistic physical systems; see the shaded intervals of 
time in Fig. 4 of the main Article. 
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